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ABSTRACT 



For processes involving structure functions and/or fragmentation functions, arguments that, over a range 
of a proper kinematic variable, there is a part that dominates the next-to-leading order (NLO) corrections 
are briefly reviewed. The arguments are tested against more recent NLO and in particular complete next- 
to-next-to-leading order (NNLO) calculations. A critical examination of when these arguments may not 
be useful is also presented. 
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1. INTRODUCTION 



In Perturbative QCD there is now a great effort towards calculating NNLO corrections 
[|l}£|]. One reason is that in several cases the NLO corrections are found to be large. Other 
reasons are that NNLO corrections are expected to increase the stability of predicted cross 
sections against changes of schemes and scales and that they will lead to more precise 
determinations of backgrounds towards uncovering signals for new physics. 

Although there is no substitute for a complete NNLO calculation, since such calculations 
are in general very complicated, as a first step one may try approximate ones. Such a step 
has been presented in || 

Below we briefly review the arguments of J3J. Sect. 2 mentions the results of certain 
more recent NLO calculations. Sect. 3 examines applications to the presently existing 
complete NNLO calculations. Finally, Sect. 4, apart from certain other points, discusses 
when approximate results may not be useful. 

For processes involving structure functions and/or fragmentation functions, in Q it was 
argued that, over a range of a proper kinematic variable, there is a part that dominates 
the NLO; and this was used to explain the fact that, in a number of the then existing NLO 
calculations, plotted against this kinematic variable, in a wide range, the cross section was 
almost a constant multiple of the Born. 

To briefly review the essential ideas of ||, consider the NLO contribution of the sub- 
prosses a(pi) + b(p 2 ) — > 7(g) + d to the large-pr process A + B — > 7 + X: 



da a 8 (ji) ^ f dx a dx b 
E-— = > / F a / A (x a ,M)F b/B (x b ,M) 

d 6 V 7T ~7 J X n Xh 



a B 6(l + + 



d 3 p 7T ah X <* X b 

+ (1-<U) (A ^B,r)^-rj), (1.1) 
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where F a /A, F h / B are parton momentum distributions to the hadrons A,B, \i and M are 
the renormalization and factorization scales, r\ the cm. pseudorapidity, 

s = (pj +p 2 ) 2 , t=(q-p 1 ) 2 , u =(q-p 2 ) 2 

and o~b and / are functions of s, t, u corresponding to the Born and the higher order 
correction (HOC). Introducing the dimensionless variables 

v = l + t/s, w = -u/(s + t) (1.2) 

(s + t + u = sv(l — w)), the HOC have the following overall structure: 

f(v, w) = f s (v, w) + f h (v, w), 



where 



f 8 (v,w) = a 1 (v)5(l-w) + b 1 (y) T — +Ci(v) 

1 - w 4- 




(a 2 (v)6(l -w) + Uv) (1 ^ w) J In ^ + c 2 (v)6(l - u;) In ± (1.3) 
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where 1/(1 — w) + and (ln(l — w)/(l — w)) + are well known distributions. The function 
fh(v,w) contains no distributions and, in general, has a complicated analytic form. 

Now denote by a s and ah the contributions of f s and fh to Eda/d 3 p and consider the 
ratio 

L = a h /(a s + a h )] (1.4) 

then, at sufficiently large Xt, for fixed total cm. energy y/S, as pr (or xt = 2pT/\fS) 
increases, \L\ decreases. 

To see the reason, consider a plot of xj, vs x a for n = (Fig. 1). The integration 
region in ( |1.1|) is bounded by w — 1, x a — 1 and xj, = 1 (hatched region). Now, for x 
not too small, F a /^(x,M) behaves like (1 — x) n ; with A =proton, n is fairly large (> 3); 
also due to scale violations, n increases as pt increases. Then contributions arising from 
the region away from w — 1 are supressed by powers of 1 — x a and/or 1 — x^. Now, in 
f s , the terms proportional to 5(1 — w) contribute at w — 1 (and so does &b) whereas 
the rest give a contribution increasing as w — > 1. On the other hand, the multitude of 
terms of fh contribute more or less uniformly in the integration region 9(1 — w) and their 
contribution ah is suppressed. As xt increases at fixed S, the integration region shrinks 
towards x a = x^ = 1 and the suppression of ah increases. 

The mechanism is tested by writing the distributions in the form @(a)] 

F a/A (x, M) = F b/B (x, M) = (1 — (1.5) 

and choosing a fictitious N » n or choosing < N << n. Then the ratio L with the first 
choice decreases faster and with the second choice decreases slower then for N = n. 

Next we neglect fh(v,w) and make the rough approximations 1/(1 — w) + ~ 5(1 — w), 
(ln(l — w)/(l — w))+ ~ 5(1 — w). Furthermore, we note that &i(t>), c\(v), 02(f), b 2 (v), 02(f) 
and part of a\(v) are either proportional to the Born term or contain the Born term times 
a smooth function of v ; the rest of a±(v) is also a smooth function of v (see e.g. Eq. (C.8) 
of @(a)] or Eq. (4.11) of @(b)]). The Born term itself is a smooth function of v. Thus as a 
first approximation we write 

f(v,w)^ Aa B (v)5(l-w) (1.6) 

where A w const. This results in Eda/d 3 p of roughly the same shape as Eda^ oxn / d 3 p 

The same argument can be made in terms of the moments of the functions 5(1 — 
w), 1/(1 — w) + , (ln(l — w)/(l — w)) + and of the functions in f h (v,w) 0(a)]. Clearly, a s 
contains all the soft, collinear and virtual contributions to Eda/d 3 p. 

At NLO the Bremsstrahlung (Brems) contributions to f s are determined via simple 
formulae Q: E.g. for gq — > 75 the Brems contributions arise from products of two graphs 
gq — > 7551. If in both graphs the emitted g arises from initial partons (g or q), the contri- 
bution in d = 4 — 2e dimensions is 

d °"'< T^,e)N c (- 2 -)(^Y\l- a r^(l + c4), (1.7) 



dvdw ' ' V £ 7 V 1 — v 
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where Tq (v,e) is essentially the Born cross section in d dimensions. If in at least one of 
the graphs the emitted g arises from the final parton (q), then 



^~7^( B , e)ClW -.(l-„)— P„(e) 



r(l-2e) f 1 



da 

where 

P«(e) = ^ J q y- £ (l - y)- £ P qq (y, e) (1.9) 

and P qq (y,e) = 2/(1 — y) — 1 — y — e(l — y), the split function in n dimensions (y < 1). 
Expanding 

(1 - „)- = - „) + ^ - e p£^>) + + (1,0) 

as well as (t> /(l — v ))~ e and t> ~ e in powers of e we determine the contributions. The singular 
terms ~ 1/e 2 and 1/e cancel by adding the loop contributions and proper counterterms. 



2. FURTHER NLO CALCULATIONS 

In addition to the examples presented in Refs. ||, the following are some NLO studies 
supporting the ideas of Sect. 1: 

(a) Large pt W and Z production in pp collisions [[J. At v^S* = 0.63 and 1.8 TeV, for 
Pt > 80 GeV the cross sections do jdp\ are also almost a constant multiple of the LO 
(Figs. 7 and 8 of §). 

(b) The production of two isolated direct photons in pp collisions [§. At VS = 1.8 TeV, 
when the pt of each photon exceeds 10 GeV the shape of the NLO QCD cross section 
da/dpT differs little from that of the Born (Fig. 2 of 0). 

Regarding NLO results for polarized reactions we mention the following: 

(a) Polarized deep inelastic Compton scattering [0, in particular the contribution of the 
subprocess jq jq to large Pt direct photon production in polarized 7 — p collisions 
(7P -> 7 + X). At y/S = 27 and 170 GeV, for x T > 0.15, it is \L\ < 0.28 and for 
sufficiently large xt, L decreases as %t —* 1 (Fig. 4 of Ref. 171]). Also, denoting by 

the 0(a>s), k = 0, 1, contributions of 7<f — > 75 to Eda/d 3 p, for 0.2 < x T < 0.8 
the factor K iq = (a^ + a ( - 1 ^)/a^ is found to differ little from a constant. 

(b) Large pr direct 7 production in longitudinally polarized hadron collisions |§||. Here 
of interest are the 0(a*), k = 1,2, contributions of the subprocess gq — > jq. As 
xt increases, the ratio —Ohjo s steadily decreases (Fig. 10 of |§). The factor K gq = 
(a^ + a^)/a^ is not constant, but increases moderately (Fig. 2 of M). 
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(c) Lepton pair production by transversely polarized hadrons fllQ|,|TT||. At fixed S, with 



increasing y/r = Mi-i+fyS, the ratio ah/cr s is found again to decrease (Fig. 3 of JT0|1 ) . 
Again, the .fT-factor is not constant, but increases moderately (Fig. 1 of ||1Q|| ). 

The considerations of Sect. 1 explain also the following fact: Taking as example large 
Pt pp — > 7 + X, at NLO, apart from the HOC of the dominant subprocess gq — > jq, there 
are contributions from the extra subprocesses gg — > qqj, qq — > qqj and qq' — > qq'j, where 
q, q' are either of different quark flavor or of the same flavor but interacting via exchange 
of a gluon. In general, these are found to be relatively small (Figs. 3, 4 and 5 of ||). The 
reason is that the extra subprocesses possess no terms involving distributions (no loops and 
vanishing contributions of the type ( |1.7| ) and (|1.8|)). 



3. NNLO CALCULATIONS 

NNLO calculations have been carried for Drell-Yan (DY) production of lepton pairs, 
W ± and Z, and for the deep inelastic (DIS) structure functions Fj(x,Q 2 ), j = 1,2 and 
the longitudinal part. Now the parts involving distributions contain also terms of the type 
(ln J (l — w)/(l — w)) + , with i = 2 and 3 and w a proper dimensionless variable. The 



subsequent calculations are carried using the updated MS CTEQ5M1 set of [fL2|, one of the 



most recent sets of NLO parton distributions We present results for /i = M = \J~Cp . 

Beginning with DY, we are interested in the process pp — > 7* + X — > + X and to 
the cross section 

da(r,S)/dQ 2 = ct(t,S) (3.1) 

where r = Q 2 / S with \/S the total cm. energy of the initial hadrons and \/Q 2 the 7* mass 
UTiJlSl. Here we deal with the subprocess q + q — > 7* and its NLO and NNLO corrections 
[ I|J] . For DY, w ~ r. We use number of flavors rif = 4. 

Denote by a^(r,S), k = 0,1,2, the 0(a k s ) part of a(r,S), by of) the part of 
arising from distributions and by the rest. Defining 

L^(r,S) = ai k \r,S)/^ k \T,S) (3.2) 

Fig. 2 shows L (fc) , k = 1,2, as functions of r for = 20 GeV. Clearly, for r > 0.3: 
L« < 0.17 and < 0.33. 

It is of interest also to see the percentage of of the total cross section determined 
up to 0(a k s ). Fig. 2 also shows the ratios 4 1) /(^ (0) + ^ (1) ) and a ( ^/(a^ + + a^) for 
the same yS; clearly, for r > 0.2 both ratios are less than 0.1. 



Now we turn to DIS [IbJTTI. Here we deal with the sum 

S(x,Q 2 ) =u v {x,Q 2 ) +d v {x,Q 2 ), (3.3) 

where u v and d v are the w-valence and d-valence quark distributions in the proton. We 
will deal with the subprocess q + 7* — > q and the nonsinglet part of its NLO and NNLO 



corrections Dill. For DIS, w ~ x. 
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Denote by ^ k \x,Q 2 ), k = 0,1,2, the 0(a k s ) contribution, by £( fe ) the part of £^ 
arising from distributions and by the rest. Defining 

LW(x,Q 2 ) = 4 fe) (x,Q 2 )/sW(x,Q 2 ) (3.4) 

Fig. 3 presents L^(x, Q 2 ), k = 1, 2, as functions of x for = 5 GeV. Now, for x < 0.5 
Z/ 1 ) is not small, but this is due to the fact that £W changes sign and £^ stays > 0, so 
at x ~ 0.3, £W vanishes. On the other hand, at x > 0.6, is less than 0.28. 

Fig. 3 also shows the ratios E?V(S (0) + £ (1) ) and £l 2) /(£ (0) + £« + £( 2 )) for the same 
Q 2 ; for x > 0.3 both ratios are less than 0.1. 

The effect of neglecting in DY or £^ in DIS is shown in Fig. 4. In DY, denoting 

K s = {aW + aP+aMya® 

K = (a® + <T W (3.5) 

we show K S (K) by solid (dashed) line at y/S = 20 GeV (upper part). Clearly, as r — > 1, 
fC s — ► i^, and for r > 0.3 the error is less than 14%. In DIS, denoting 

K s = (£<°> + )/E(°) 

K = (E(°) + EW+S (2) )/S (0) , (3.6) 

we show and K at v/Q 7 = 5 GeV (lower part). Again, as x — > 1, K s — > K. Now, in 
spite of the fact that is, in general, not small, K s differs from K even less. The reason 
is that the NLO and NNLO corrections are smaller than in DY, and so are E^ /£*-°). 

4. CONCLUDING REMARKS 

The above discussion and examples show that for processes involving structure and/or 
fragmentation functions, for not too small values of a proper kinematic variable (xt for 
large-pr reactions, r for DY, x for DIS), one may retain only that part of the differential 
cross section arising from distributions [jlg|| . At NNLO the range of this variable is larger 
for DY than for DIS (Figs. 2 and 3). Yet, regarding the .fT-factor, which determines the 
physically important quantity, DIS is somewhat advantageous (Fig. 4). 

As we go to NNLO, in view of the presence of terms of the type (ln l (l — w)/(l — w)) + 
with i = 2 and 3, the rough approximation of replacing (ln(l — w)/(l — w))+ by 5(1 — w) 
is becoming worse. In general, this implies that with NNLO corrections, the shape of a 
physical quantity should deviate more than that of the Born term. 

It is desirable (and nontrivial) to extend the formulas (|1.7|) and (|1.8|) to the NNLO and 
perhaps even higher orders. 

The question now is when the arguments of Sec. 1 may not be useful. Such a case 
is when, over a wide range of w, is comparable and of opposite sign to cx( fe) . Then 
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r ( fc ) _|_ ff w i s small and is large in absolute value. Even then, for w very near 1, 



should decrease, but in that case threshold resummation |19],^0| is important, and the 
approximation is not useful. Of course, in such a case, the correction \a^\ = + cr^| 
will be small. The point, however, is that we do not see how one can determine such a case 
without calculating a h . 
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FIGURE CAPTIONS 



Fig. 1 The integration region in the expression ( |1 . 1| ) for cm. pseudorapidity rj = 0. 

Fig. 2 The ratios = a® /<r<$ and of '/{a^ + <x« + a^) (solid lines) as well as L« = 
(7^ /(7^ and cjj, / (cr^ +a^) (dash-dotted lines) for Drell-Yan lepton-pair production 
versus r = Q 2 /S at VS = 20 GeV. 

Fig. 3 The ratios = sf } /S (2) and E ( f /(E^+E^+E^) (solid lines), where £ = u v +d v , 
as well as L« = Ej^/EW and E [ ^/(E^ + £ (1) ) (dash-dotted lines) for q + 7* -> g 
versus 2 at Q 2 = 25 GeV 2 . 

Fig. 4 .fT-factors: Approximate K s (solid lines) and exact K (dotted lines) for the Drell-Yan 
case of Fig. 2 (upper part) and for the case of Fig. 3 (lower part). 
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Fig. 3 




Fig. 4 



